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ABSTRACT 

We present a model for high-energy emission sources generated by a standing magnetohydro- 
dynamical (MHD) shock in a black hole magnetosphere. The black hole magnetosphere would 
be constructed around a black hole with an accretion disk, where a global magnetic field could 
be originated by currents in the accretion disk and its corona. Such a black hole magnetosphere 
may be considered as a model for the central engine of active galactic nuclei, some compact X-ray 
sources and gamma-ray bursts. The energy sources of the emission from the magnetosphere are 
the gravitational and electromagnetic energies of magnetized accreting matters and the rotational 
energy of a rotating black hole. When the MHD shock generates in MHD accretion flows onto 
the black hole, the plasma's kinetic energy and hole's rotational energy can convert to radiative 
energy. In this letter, we demonstrate the huge energy output at the shock front by showing 
negative energy postshock accreting MHD flows for a rapidly rotating black hole. This means 
that the extracted energy from the black hole can convert to the radiative energy at the MHD 
shock front. When axisymmetric shock front is formed, we expect a ring-shaped region with very 
hot plasma near the black hole; the look would be like an "aurora" . The high energy radiation 
generated from there would carry to us the information for the curved spacetime due to the 
strong gravity. 

Subject headings: accretion, accretion disks — black hole physics — magnetohydrodynamics 
(MHD) — shock waves 



Introduction 



A number of observations provides evidence for super-massive black holes in galactic centers or stellar- 
mass black hole in galaxies (e.g. Mivoshi et al. 1995 : Ghez et al. 1998 . 2005t Genzel et ahlbood : Schodel et al 



2003: Bender et al. 2005: Eisenhauer et al 



20051 ). Several observations show high-energy phenomena around 



black hole candidates where some amount of energy in accretion flows is released. In this paper, to understand 
the high-energy phenomena near the black hole and the nature of the curved spacetime due to the hole's 
strong gravity, we consider ingoing magnetohydrodynamic (MHD) flows in the magnetosphere of a black 
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hole, which is considered as the transit region from an accretion disk to the event horizon. The global 
magnetic field in the magnetosphere could be generated by currents streaming in the disk and its corona 
around the black hole. 

The s tructure of a steady and a xisym metric black hole magnetosphe re with a thin di sk has been dis- 
cussed by iTomimatsu fc Takahashil (120011) as a vacuum magnetosphere . lUzdenskvl (|2005l ) as a force- free 
magnetosphere, and ICamenzindl (|l987l) : iNitta. Takahashi. fc TomimatsrJ (|l99l[ ) as a MHD magnetosphere. 
These models show the disk-black hole connecting magnetic field lines, as shown in Figure [T^. The loop-like 
shaped magnetic field is generated around the disk's inner edge region. Some magnetic field lines connect 
to the event horizon. Due to the plasma inertia effect, the magnetic field lines are forced to direct toward 
the event horizon, and then almost radial flows would be reasonable at least near the event horizon. Such 
a configuration of the loop-like magnetic field lines indicates the interaction between the black hole and the 
disk surface (or corona) located at several times the inner-edge radius of the thin disk. Along this loop- 
like magnetic field lines, the MHD fiuid ejected from the disk surface streams inward due to the dominant 
gravitational force by the black hole. 

Although in many previous works equatorial accretion flows have been discussed, by considering the 
black hole magnetosphere, we expect naturally off-equatorial inflows streaming along the lo op-like magnetic 



field lines. So, we also expect the shock formation near the polar region of the black hole fsee lFukumura et al 



20071 for adiabatic shocks). Of course, the shock formation by equatorial accretion flows would be possible 
within the required physical conditions. However, the emission from the off-equatorial MHD shock would 
include the information for the "magnetosphere" (i.e., the configuration of the magnetic fields) around a 
black hole, and would enable us to find the information for the black hole spacetime related to the radiation. 

We are interested in the energy conversion from the accreting plasma's energy, which is composed 
of the kinetic part and magnetic part, to the radiative energy at the inner-most region of the black hole 
magnetosphere. This is because such energy conversion by the shock formation would be related to the 
generation of huge radiative energy in the hole's deep gravitational potential. Furthermore, in addition to 
the plasma's kinetic energy, the rotational energy of the black hole may also be converted to the radiative 
energy (see Fig.[TjD). For a fast MHD shock on accretion onto a black hole, about 10% of the upstream flow's 
energ y (including the rest m ass energy of a particle) can convert to the thermal energy of the downstream 



flow ( Takahashi et al.l 20061 ). This means that the plasma temperature becomes so high across the shock 



front. Such a hot plasma region can be considered as a source of high energy radiation, which gives us both 
the information for the strong gravitational field and the state of magnetized plasma around the black hole. 
Of course, some part of the radiation emitted from the hot plasma will fall into the black hole because of the 
strong gravitational lens effects. Nevertheless we can expect that huge radiative energy can be released at 
the very hot shocked plasma region, and considerable radiation flux will be obtained for us. 

In this paper, we discuss energy release by shocked MHD accretion onto a rapidly rotating black hole. 
In §2, we introduce trans -fast magnetosonic accretion flow s and apply the shock condition to the flows. Here, 
we adopt the method by Takahashi fc Tomimatsu ( 20081 ) for solving trans- fast magnetosonic flow solutions. 
Then, we can easily obtain the solutions without critical condition analysis at the magnetosonic points. In 
§3, we show negative energy postshock MHD inflows that means the energy extraction from a rotating black 
hole and discuss the released energy at the shocked region. In §4, we summarize our MHD shock models. 



^ Th e escape rate of photons from the ergosphere is almost 10% and this rate becomes larger for a larger spin value (see 
iTakahashi fc Takahashi.20ia) . 
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Fig. 1. — (a) A model of the black hole magnetosphere. Ingoing MHD flows stream along the disk-black hole 
connecting magnetic field lines. Because of the loop-shaped magnetic field configuration, off-equatorial shock 
would be expected, (b) MHD shock in black hole accretion. After passing through first fast magnetosonic 
point, the MHD inflow can make the MHD shock. The postshocked sub-fast magnetosonic flows must pass 
through the second fast magnetosonic point again before reaching the event horizon. 



MHD Accretion with MHD Shock 



We consider MHD flows in a stationary and axisymmetric magnetosphere around a rotating black hole. 
The background metric g^^i, is written by the Boyer-Lindquist coordinate with c = G = 1. The basic 
equations for MHD flows are the equation of the particle number conservation, the equation of motion and 
Maxwell equations. We also assume the ideal MHD condition and the polytropic relation for the plasma 
flows. There are five field-aligned flow parameters on the fiows; that is, the total energy E, the total angular 
momentum L, the angular velocity o f t he magnetic fleld line ilp and the number flux per flux-tube rj and 
the entropy S (see Camenzindl 1986al lb: Takahashi et al. IQQOl for the definitions of these flow parameters). 
The accreting flow onto a black hole must pass through the slow magnetosonic point, the Alfven point, and 
the fast magnetosonic point, in this order. The critical condit ions at the Alfve n point and the magnetosonic 
points restrict the acceptable ranges of these parameters (see lTakahashill2002L for the details). The poloidal 



r.2 



(Fe0)2]/p^, and the 
-I- cos^ 9, 



magnetic field Bp seen by a lab-frame observer is defined by Bp 

toroidal magnetic field is defined by B^ = {A/Y,)Fre, where A = r'^ — 2mr + a'^, S 

— A sin^ 9, and m and a are the hole's mass and spin, respectively. The term i^^i^ is the electromagnetic 
tensor. In the following, we assume the cold inflows (S — 0) except for the shock front. Although, at the 
shock front, the plasma temperature can become so high, it falls off immediately by the radiation. 

The relativistic Bernoulli equation for cold MHD flows, whi ch determines the poloidal velo city (or the 
Alfven Mach number) along a magnetic tube, can be written as ( Takahashi fc Tomimatsu 20081) 



-a- M\aBl+Bl) ^0 , 



(1) 



where e = E ~ QpL with E = E/fXc and L = L/fic- The enthalpy for cold fiows is denoted by fXc = Wpart, 
where mpart is the particle's mass. The terms Bp = Bp/{4:T:ficV) ^'t> = B^/ [A-k^cIIPw) are introduced 
to non-dimensionalize, and the latter is given in terms of the relativistic Alfven Mach number and the flow 
parameters; 

GtL 



Bi 



G^E 



(2) 
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where Gt 



9tt 



M is defined by AP = {ul/u^ 



defined as 



'■AW 



Takahashi et al 



gt^rip, Gel, = gttj, + gijicjiilF and a = Gt + G^flp- The relativistic Alfven Mach number 
P = (Wp/wAw)'^ ~ '^'p/^p^ wliere Up is tlie poloidal velocity, the A lfven velocity haw is 
aBp/{ATT^cn), and n is the number density of the plasma (see also Camenzind 1986al lbl: 



1990l) 



Differentiating eqs. ([TJ and ([2]) along a magnetic field line, we obtain the velocity gradient at any point, 
and see the critical point in flow solutions. That is, the flow solution becomes singular at the magnetosonic 
points and the Alfven points without the regularity conditions. To obtain a physical solution tha t is smooth 
every where, the complicated parameter search for satisfying regularity conditions is required (jTakahashi 
20021 ). To avoid t he heavy task, we will introduc e the ratio of the poloidal and toroidal components of the 
magnetic field as (jTakahashi fc TomimatsTjl2008l ) 



(3) 



to specify a magnetic flux-tube for the purpose of solving the relativistic Bernoulli equation ([T]). That is, 
the poloidal magnetic field is related to the toroidal magnetic field through the function (3{r,d), which is 
assumed to be regular at magnetosonic points, where we do not assume the concrete function for the poloidal 
magnetic field. We eliminate Bp in eq. ([T]) by eq. and also eliminate B"^ by eq. ([2]). Then, eq. ^ can be 
reduced to the quadratic equation of AP. As a result, we obtain trans- magnetosonic solutions easily without 
the analysis for the critical conditions of the magnetosonic points. Although the toroidal magnetic field has 
several constraints at some points along the flow (i.e., the event hor izon, the A l fven p oint, the Anchor point, 
etc), the typical feature of for a given Bp has been discussed in iTakahashil (|2002l ) by using the standard 
critical point analysis. So, we can apply it to make a reasonable test function for B^{r, 6). For example, we 
can introduce it as a simple form of 



= i-g^^mF - UJH? [l + C (A/E) 



(4) 



where C is a constant. The parameter C is introduced to specify the discontinuity of _B|| at the shock front, 
where B\\ is the components of magnetic field parallel to the shock front. 

Now, we will consider shocked MHD accretion onto a black hole. To obtain a shocked black hole accretion 
solution, it is necessary to set up two trans-magnetosonic solutions, which correspond to the upstream and 
downstream solutions. At the shock location, where the shock conditions are required, the upstream super- 
fast magnetosonic solution is connected to the downstream sub-fast magnetosonic solution within suitable 
ranges for the field-aligned flow parameters. After the shock, the sub-fast magnetosonic inflow passes through 
the second fast magnetosonic point, and then it falls into the horizon. In general, the downstream solution 
would have different values for the field-aligned parameters from the upstream solution due to the radiation 
loss at the shock front. Hereafter, we assume the conservation of the angular momentum of the flow L, the 
angular velocity Vlp oi the magnetic field line and the number flux r] across the shock front, while the energy 
of the fiow jumps across the shock front; that is E^^p > -Edown- 

The conditions for the MH D shock formation in the general relativistic framework are summarize by. 



Lichnerowicz ( 1967 . 1976[) : that is, t he particle number, the energy momentum and the mag netic flux 



e-g-, 

are conserved across the shock (see also lAppl fc Camenzindl 1 19881 : iTakahashi et al.l l2002l |2006[ ). For an 
isothermal MHD shock in cold plasma flow, from the jump condition for the radial component of the energy 
momentum conservation, we obtain the simple relation of 



1 



AP 



1 



(5) 
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Fig. 2. — The Alfven Mach number for shocked MHD accretion onto a rotating black hole with a — 0.8m, 
r^F = 0.4rii„ax, LVlp = —0.5 and 9 = tt/2. The preshock solution of Eup = 1-0 is plotted by the red thick 
curve, while the postshock solutions for i?down = 0.8 (S) and i?down = —0.43 (S') are shown by the blue 
curves. As a parameter of the magnetic field line, we set Cup = —5.0 and Cdown = 40.0. The broken curves 
show the Alfven Mach number related to the fast magnetosonic wave speed for the upstream and downstream 
flows. The dotted line labeled sk shows the Alfven Much number M^^^^ given by eq.®. The shock radius 
for each postshock solution is shown by the intersection between a blue curve and the dotted black line. 



where for the sake of simplicity we assume that B^^ 



Ddown 



\B' 



up 



< 



I Ddown 



(for the fast-magnetosonic 

shock) and the shock front is normal to the upstream flow in the poloidal plane. The value for 0^ (i.e., the 
value for C) increases across the fast magnetosonic shock, and then the Mach number decrease; that is, the 
kinetic energy of the flow is released there. 



High Energy Radiation powered by Rotating Black Hole 



Rotational energy of a rapidly rotating black hole can be extracted by MHD inflows wi th negative energy 
due to the interaction betwe en magnetic and gravitational fields in the black hole spacetime (jBlandford fc Znaiek 



19771: iTakahashi et al.lll990D . The extracted energy is carried to the magnetosphere in the form of outgoing 
Poynting flux. If the extracted energy is converted to some kinds of fluid's energy related to the radiative 
process directly, it may be observable for us. However, in the framework of original Blandford-Znajek (BZ) 



- 6 - 




Fig. 3. — The relation between the energy of the postshock flow and the shock location for the solutions for 
shocked MHD accretion onto a rotating black hole with a — 0.8m and Vlp = 0.4ilinax (thick curves). The 
values of Cdown are selected from 0.0 to 100.0, while Cup = —5.0. The shaded region is the forbidden region 
for the fast MHD shock. The broken curve (Adown) shows the location of the Alfven point for the postshock 
flow, whose location depends on the energy i?down ■ The vertical line shows the location of fast magnetosonic 
point Fup of the preshock flow. 



process in the force-free magnetosphere, realistic conversion mechanisms from the magnetic energy to the 
other form are still not clear, because in the force-free approximation the inertia of matters is ignored. 

We are now discussing MHD accretion onto a black hole. The MHD inflow takes the fluid energy into 
the hole (that is, the kinetic energy, the internal energy and the rest mass energy of the fluid), which is 
positive at the plasma source. However, the situation of energy extraction from a rotating black hole can be 
also described in the MHD magnetosphere. In this case, the negative energy MHD inflow {E < 0) is realized, 
where it is required that the Alfven point of the flow considered is located close to the inner light surface in 
addition to the conditi on for the BZ process; both the Alfven point and the inner light surface need to posit 
inside the ergosphere (jTakahashi et al.l 119901 ). By considering negative energy MHD inflows, the rotational 
energy of the black hole can be carried in the surrounding magnetosp here as the outward energy flux to the 
distant region directly (^McKin nev fc Gammiel 2004 ; McKinnev 20061) or indirectly by way of an equatorial 
disk (iLiii2002c IWang et al..,2003.) . In the latter case, the extracted energy from the hole is deposited onto the 
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disk, and then it converts to the thermal energy of the disk plasma and/or the kinetic energy of the outgoing 
disk winds. 

Now, at the fast MHD shock, we understand that the kinetic energy of the upstream MHD flow converts 
to the thermal and magnetic energies of the downstream flow. The most efflcient energy conversion to 
radiation would be achieved by isothermal shock; where it is assumed that the thermal energy generated at 
the shock front escapes from the shocked plasma immediately. When the negative energy accretes onto the 
black hole as the postshock flow, we can also expect that the extracted hole's energy would convert to the 
radiative energy at the shock front (see Fig.[T}D). By considering the regularity condition at the Alfven point, 
which is related to the amount of the jump of energy and angular momentum be tween the preshock and 



posts hock solutions, we find the necessary condition for such a energy release process (jTakahashi fc Takahashi 



201(1 in preparation). 

Figure [5] shows two examples of solutions for shocked accretion onto a black hole. After passing through 
the Alfven point Aup and fast magnetosonic point Fup, the accretion flow becomes super- fast magnetosonic, 
and then it falls into the black hole. When the MHD shock arises at the radius labeled S or S' on the 
accreting flow (see the vertical downward-arrows), the postshock flow becomes sub-fast magnetosonic. Then, 
the postshock flow must become super-fast magnetosonic again by passing through the fast magnetosonic 
point Fdown or F^j^^^^ located between the shock front and the event horizon. Note that the branch of 
the postshock solution labeled S connects to the Alfven point Adown, and it is also a shock- free accretion 
solution. On the other hand, the postshock flow solution labeled S' shows an unphysical solution as a shock- 
free solution, because the branch of the solution does not connect to the Alfven point and diverges outside 
the shock location. We can not accept such a divergent solution as a shock-free solution streaming from the 
plasma source to the horizon. However, by considering the shock formation, this solution S' can survive as 
a physical postshock flow solution in the region within the shock radius. 

At last, we consider the conditions for the formation of the postshock solution with negative energy and 
the possible range of the amount of the energy extracted from the black hole. Figure |3] shows the relation 
between the energy of the postshock flow E^own and the shock radius for the accretion flow onto the 
black hole with a = 0.8m. The energy of the preshock flow solution is set to E^p = 1.0 and the angular 
velocity of the magnetic field lies is set to flp = 0.4r2„iax, where Omax is the maximum value for the field line 



exist ing the plasma source region located between the inner and outer light surfaces (see I Takahashi et al 



19901 ) . Some thick curves in Figure |3] correspond to the different values of Cdown for the postshock solution, 
and we set Cup = —5.0 for the preshock flow solution. Although we consider the case where the shock 
front is generated inside the ergosphere, the postshock flow with negative energy {E < 0) is possible. This 
means the extraction of the hole's rotational energy by the ingoing MHD flow. The minimum energy of the 
postshock flow i?niin IS given by E = LVlp < 0. The extracted energy can be released at the shock front 
where the hot plasma region is generated. From the hot plasma in the shock front, the observable energy 
output is expected; e.g., in the form of high-energy radiation, bulk- motion of outgoing plasma flows, etc. 



4. Concluding Remarks 

We have discussed global shocked accretion solutions in a black hole magnetosphere, which are composed 
of two trans-magnetosonic solutions with a MHD shock where the general relativistic MHD has been applied 
to the accretion flows. Then, we have shown the possibility of a very hot shocked plasma region very close 
to the event horizon. At this stage, in the study of the black hole magnetosphere, the expected energy 
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spectrum from the MHD shock is not clear. Nevertheless, we newly point out the possibility of the energy 
extraction from the black hole through the MHD shock with the negative energy in the postshock flow. 
Although, in this paper, we show the possibility of shocked accretion solutions with negative energy {E < 0), 
the dependence of flow parameters (i.e., -Eup, ^f, o,, ^up) and flow's latitude will be presented in the next 
paper (Takahashi & Takahashi, in preparation). 

To complete the black hole's rotational energy release at the MHD shock, we must investigate the 
structure of the shocked hot plasma region. Although the value of the released energy depends on the 
field-aligned fiow parameters, we can estimate the plasma density and temperature there by the help of the 
accretion disk models that give the boundary condition for the magnetosphere. We need a realistic model of 
the black hole magnetosphere with a magnetized accretion disk. Although we assume a given magnetic field 
line, the configuration may be unstable. Then, the dynam ical behaviors of the field configuratio n should 



be also investiga ted by helping numerical simulations (e.g.. iKoide et al.ll2006t iMcKinnev fc Gamm ic 2004j; 



McKinnevI 120061 ) . Then, a significant amount of X-rays and gamma-rays would be estimated. In order to 
take into account such radiation effects in the general relativistic plasma, we should formulate the local 
radiation effects of MHD fluid by using the energy-momentum tensor with the local heat flux, in addition 
to the fluid part and the electromagnetic part. Then, we will find the evidences for real black holes in the 
observational data (e.g., the spectrum in high-energy range, the short time variability, the radio image of 
black hole shadow. X-ray polarization, etc). That is, we expect that the high energy emission from this 
hot plasma brings to us additional information about the black hole spacetime; especially, the polar region 
emission including this information can be distinct from the emissions from the equatorial plasma source, 
which have been investigated by many authors in the models of accreting gas disk. 

We are grateful to Sachiko Tsuruta for her helpful comments. This work was supported in part by the 
Grants-in-Aid of the Ministry of Education, Culture, Sports, Science and Technology of Japan [19540282 
(M.T.), 21740149 (R.T.)]. 
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